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ABSTRACT 

B2013-I-370 and B2023+336 are two blazars at low-galactic latitude that were previously proposed 
to be the counterparts for the EGRET unidentified sources 3EG J2016-I-3657 and 3EG J2027-h3429. 
Gamma-ray emission associated with the EGRET sources has been detected by the Fermi Gamma- 
ray Space Telescope, and the two sources, IFGL J2015. 7-1-3708 and IFGL J2027.6-(-3335, have been 
classified as unidentified in the 1-year catalog. This analysis of the Fermi-LAT data collected during 31 
months reveals that the IFGL sources are spatially compatible with the blazars, and are significantly 
variable, supporting the hypothesis of extragalactic origin for the gamma-ray emission. The gamma- 
ray light curves are compared with 15 GHz radio light curves from the 40-m telescope at the Owens 
Valley Radio Observatory (OVRO). Simultaneous variability is seen in both bands for the two blazar 
candidates. The study is completed with the X-ray analysis of IFGL J2015. 7-1-3708 using Swift 
observations that were triggered in August 2010 by a Fermi-detected flare. The resulting spectral 
energy distribution shows a two-component structure typical of blazars. We also identify a second 
source in the field of view of IFGL J2027. 6-1-3335 with similar characteristics to the known LAT 
pulsars. This study gives solid evidence favoring blazar counterparts for these two unidentified EGRET 
and Fermi sources, supporting the hypothesis that a number of unidentified gamma-ray sources at 
low galactic latitudes are indeed of extragalactic origin. 

Subject headings: Galaxies: individual (B2013-I-370, B2023-I-336) - Gamma rays: galaxies 



1. INTRODUCTION 

Blazars are a subclass of Active Galactic Nuclei (AGN) 
which host supermassive black holes in their centers and 
have relativistic jets pointing along the line of sight to 
the observer. They show broadband non-thermal emis- 
sion extending from radio frequencies to gamma rays, 
with a typical two-component spectral energy distribu- 
tion (SEE), see, e.g. 
Fermi Larg e Area Te, 



are detected at high G a.lactic latitudes, and only 14 are 



Sambruna et al. 1996). In the 
es cope [F ermi-LAL') 1-year cata- 
log (IFGL, |Abdo et al.||2010b[) 573 out of 1451 (39%) 



detected objects are associated with known blazars, mak- 
ing them the most numerous class of gamma-ray sources. 
Only a few bright blazars have shown gamma-ray vari- 
ability correlated with other wavelengths, allowing a firm 
identification. Almost all blazar associations reported in 
the IFGL catalog are based on positional coincidence of 
the Fermi-hAT gamma-ray source with the radio loca- 
tion of a preselected blazar. Most gamma-ray blazars 
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reported at |&| < 10° (Abdo et al. 2010b). The same 



deficit of detections at low Galactic latitudes is reported 
for gamma-ray emitting AGNs, which includes blazars 
and also misaligned a ctive galaxies. The first Fermi-hAT 
AGN catalog (lLAC, [Abdo et al.|2010c[ ) makes 51 associ- 
ations with low-latitude AGJN s from the VLBA calibrator 
survey, but a simple extrapolation assuming an isotropic 
distribution of gamma-ray emitting AGNs indicates that 
Fermi-LAT should be detecting - 150 AGNs at |6| < 10° 
( |Abdo et al.]|2010c[ ). 

'I'his spatial anisotropy in the distribution of gamma- 
ray blazars and AGNs is mainly caused by a lack of extra- 
galactic source catalogs covering low Galactic latitudes. 
Catalogs of candidate gamma-ray blazars are generally 
built by selecting compact flat-spectrum radio sources 
and identifying and classifyin g them as blazars t hrough 
optical spectroscopy (see, e.g. Healey et al.||2008 ). Close 
to the Galactic plane, diffuse radio emission and confu- 
sion with local radio sources make candidate blazars dif- 
ficult to s elect. Heavy optical extinction due to interstel- 



lar dust (Schlegel et al 
studies at low Galactic 



1998| further complicates AGN 
The three pre-selected 



atitudes 

catalogs used to find blazar associations for IFGL sources 
avoid th e Galactic plane due to these observational dif - 



20071 [20081 [Massaro et al.||2009 ) 



ficulties ( Healey et al. 

Therefore, the deficit of gamma-ray blazars at low Gaiac- 
tic latitudes is mainly an artifact of the association pro- 
cedure caused by our limited observational knowledge 
of blazars behind the Galactic plane. In fact, unbiased 
gamma-ray surveys like the Fermi 1-year catalog might 
be the best tools to identify the population of blazars 
at low Galactic latitudes that are hardly visible at other 
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TABLE 1 

NON-IFGL SOURCES INCLUDED IN THE MODEL 



Name 


RA 


DEC 


<^95% 


TS 


Fit 


r(Sc) 


Flux 1-100 GeV 


2FGL 




[dog] 


[dog] 


[deg] 








[10^9 s-i] 


association 


J2020. 0+4158 


304.70 


42.16 


0.08 


232 


PL 


2.94 ± 0.03 


3.50 ± 0.22 


J2020.0+4159 


J2022. 2+3840 


305.55 


38.68 


0.13 


117 


PL 


3.13 ± 0.04 


2.45 ± 0.16 


J2022.8+3843cFl 


J2017.4+3628 


304.45 


36.42 


0.02 


465 


PL 


2.59 ± 0.02 


7.72 ± 0.28 


J2018. 0+3626-^ 


J2025. 1+3342 


306.26 


33.70 


0.05 


289 


PL 


2.94 ± 0.03 


3.25 ± 0.14 


J2025. 1+3341 


J2028. 3+3333 


307.08 


33.55 


0.03 


1112 


PL Exp 


0.97 ± 0.05 


10.74 ± 0.38 


J2028. 3+3332 














(1370 ± 31 MeV) 







wavelengths. 

The third EGRET Catalog ( |Hartman et al.|1999[ ) listed 
66 high-confidence associations with blazars, but none 
of them is at |6| < 10°. Of the 80 reported sources 
at low Galactic latitudes, 74 were left unidentified. A 
re-evaluation of the correlation of fiat-spectrum radio 
sources with northern EGRET sources suggested blazar 
identificat ions for 5 lo w-latitude unidentified sources 
( [Sowards-Emmerd et al.|200 3 ) . Later detailed studies on 
individual EGKE'i' tinidentified sources suggested blazar 
counterparts for 3EG J2016-h3657 and 3EG J2027+3429. 
In this paper we establish firm blazar associations for 
these two unidentified gamma-ray sources. 

The blazar B2013-I-370 (RA = 20'^15™28.80^ DEC = 
37°10'58") is located at Galactic latitude b = 1.22° in the 
Cygnus region. It was proposed as a counterpart for the 
EGRET unidentified source 3EG J2016-F3657 in a com- 
prehensive study of the X-ray sources in the field followed 
by optical spectroscopy of the candidate counterparts 
( [Mukherjee et a l. 2000; Hal pern et al.|2001[ ). B2013+370 
was originally detected as a fiat-spectrurij radio source 
( |Duin et al.]|l975[ [Weiler fc Shaver] | 1978[). It appear s 
compact down to milliarcsecond scales ( |Lee et al .|2008[) 



There is no spectroscopic redshift, although the opti- 
cal counter p art is manifest via its variability (Halpern 



et al. 20011. Massaro et al. (20091 classify B2013+370 



as a b lazar of uncertain type, and |Veron-Cetty fc Veron| 



(2010) list it as a presumed BL Lac object. Very High 
Energy (VHE) gamma-ray observations by the Whipple 
10-m telescope led to a 99% confidence level (c.l.) flux 
upper limit of F{E > 350 GeV) < 2.0 x 10"" cm-^ g-i 
assuming a photon index of 2.5 ( Fegan et al.| |2005). 

B2023+336 (RA = 20''25'°10.84^ DEC = 33°43'00"), 
another blazar in the Cygnus region located at 6 = 
—2.37°, was also propo sed as counterpart for an EGRET 
unidentified source. Sowards-Emmerd et al.] (20031 
first suggested the association ot 3EG J 2027+3429 with 
B2023-I-336 after re-evaluating the correlation between 
EGRET unidentified sources and flat-spectrum radio 
sources. A later study supported the association by re- 
porting sign iflcant flux variabih ty from the blazar in the 
X-ray band (S guera et al.|200"4 ) . Previously, other stud- 
ies had suggested a young pulsar as counterpart for the 
EGRET source given its positional coincidence with an 



OB association (Romero et al.]]1999 



Zhang et al.]]2000| . 

At lower energies, B2023-I-336 has been detected as a 
compact flat-sp e ctrum radio source (j Criffith et al.|]T990 



Lee et al.|2008| ) . Sowards-Emmerd et al. ( 2003^) obtamec 



a spectrum of the optical counterpart coincident with 
the radio position. The authors classified B2023+336 as 
a BL Lac object and derived a redshift oi z — 0.219 based 



on H spectral lines. However, the signal-to-noise ratio of 
the optical spectrum is quite low due to heavy optical ex- 
tinction, and the obtained redshift should be considered 
as tentative. 

In the Fermi 1-year catalog, the blazars B2013-I-370 
and B2023+336 are respectively found in the vicin- 
ity of the unidentifie d sources IFGL J20 15. 7+3708 and 
IFGL J2027.6+3335 ( ]Abdo et al.]]2010b] ). 

In this paper we determine, discuss, and establish 
firm blazar associations for the gamma-ray sources 
3EG J2016-I-3657 and 3EG J2027+3429 based on new 
gamma-ray. X-ray, optical and radio data. 

2. GAMMA-RAY OBSERVATIONS 

The LAT is the main instrument on-board the Fermi 
Gamma-ray Space Telescope. It is a multi-purpose tele- 
scope that surveys the gamma-ray sky between 100 MeV 
and 300 GeV, including the largely unexplored energy 
window between 10-100 GeV. 

We selected photons from the first 31 months of the 
Fermi mission (Aug 2008 - Feb 2011). The analysis was 
performed using the Fermi-hAT ScienceTools software 
pack age (version v9rl8). D iffuse class photons (Pass6 V3 



IRE, [Atwood et al.]]2009| | with E > 300 MeV were se- 
lected in a 17" x 17" region around B2013+370, a Region 
of Interest (ROI) which also contains B2023-I-336. 

We generated a model that includes all the sources 
in the Fermi 1-year catalog inside the ROI (46 IFGL 
sources total), and the Galactic and isotropic diffuse 
gamma-ray emissior[^ The IFGL sources were modeled 
by a simple power-law spectrum of the form: dN/dE = 

No{E/Eo)^^ where the photon index F and the nor- 
malization iVo were left free. The only exc eption is 
for the known gamma-ray pulsars in the ROI ( |Abdo et 
al. 2010d), which have been described by a power-law 
with exponential cutoff, dN/dE = Nq {E/Eo)~^ (T^I^-^ 
where the cutoff energy, i^c, was also left free. Known 
gamma-ray sources outside the ROI but within 5° from 
the edges were included in the model, but their spec- 
tral parameters were fixed to the IFGL values. The 
isotropic diffuse background was modeled using the pub- 
lic spectral template (isotropic_iem_v02.txt) with a free 
normalizatior0 Galactic diffuse emission was also mod- 
eled by a template released by the Fermi Collaboration 
(gll_iem_v02.fit)®. Because the sources of interest lie close 
to the Galactic plane and the background is therefore 

* This analysis was completed before the Fermi 2-year catalog 
was released, and only information from IFGL was included in the 
analysis. 

^ http : //f ermi .gsfc. nasa.gov/ssc/ 
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TABLE 2 

Analysis results for the gamma-ray sources studied in this work 



Name 


RA 
[deg] 


DEC 
[deg] 


<^95% 
[deg] 


TS 


Fit 


r(Sc) 


Flux 1-100 GeV 

[10-'' cm-2 s-i] 


IFGL J2015. 7+3708 
J2025. 1+3342 
J2028.3+3333 


303.89 
306.26 
307.08 


37.17 
33.70 
33.55 


0.03 
0.05 
0.03 


808 
289 
1112 


PL 
PL 
PL Exp 


2.57 + 0.02 
2.94 + 0.03 
0.97 + 0.05 (1.37 + 0.3 GeV) 


9.26 + 0.26 
3.25 + 0.14 
10.74 + 0.38 





I >« 

f y 



3EG J2021+3716, 



1 FGL J2030.0+3641 



3EG J2027+3429 



1FGLJ202l%j9l\ 



1FGL J2015.'?v.3708 



3EG J201 6+3657 



1FGLJ2027.6+3335^ B2023+336 

°° J2025.1+3342 
J2028.3+3333 



NVSS{1.4. 




1FGLJ2027.6+3335 

oO 

J2028.3+3333 



J2025. 1+3342 



Fig. 1.— Top: Fermi-LAT counts map {E > 300 MeV) cov- 
ering B2013+370 and B2023+336 where the galactic and diffuse 
gamma-ray emission have be en subtracted. Bot tom: Closo-up 
NVSS counts maps (1.4 GHz, [Condon et al.||1998| l of B2013+370 
(left) andB2023+336 (right). White dashed con tours indicate 95% 
c.l. error ellipses for EGRET sources from Ma ttox et al.||2b01| . 
The error ellipses for the sources in Fermi 1-year catalog are shown 
in cyan. Green circles indicate the 95% c.l. error circles derived in 
this analysis, and yellow markers indicate the catalog positions of 
the two blazars. 

quite heterogeneous, we weighted the Galactic diffuse 
emission by a power-law with free normalization and in- 
dex, yielding a decrease of 11% on the normalization and 
a spectral hardening of 3% with respect to the public 
template. These additional degrees of freedom signifi- 
cantly improved the overall likelihood of the fit. 

A binned likelihood analysis was performed using the 
gtlike tool. Using a residuals map, we identified ad- 
ditional sources not in the 1-year catalog, and included 
them in the model if they had TS > 25 , equivalent to 



a significance cr ~ 5 (Mattox et al. 19961. The location 



of the additional sources was refined using gtf indsrc 
De tails on these additional sources are included in Ta- 

To derive the spectra and light curves of the studied 



sources the data were divided in time bins for the light 
curves and energy bins for the spectra, and the like- 
lihood analysis was performed on each bin. To avoid 
having too many free parameters, the Galactic diffuse 
emission was fixed. The photon indices of all sources 
were also fixed, but their normalizations left free. The 
spectra were generated using the likeSED python mod- 
ulcp^ Spectral points in the spectra were required to 
have TS > 9. To keep the maximum number of data 
points in the light curve s, flux points were d erived when 
TS* > 1 as described in |Abdo et al.| ( |2009a| ). Otherwise 
90% confidence level upper limits were calculated. 

In our analysis we find that a single source at the loca- 
tion of IFGL J2027. 6+3335 does not represent the data 
well, and so it is removed from the model and replaced 
by two separate point-like sources (see Figure fl]). A like- 
lihood ratio test (LRT) confirms that the data favors the 
2-source hypothesis. We compared a model with only 
the Fermi 1-year catalog source, IFGL J2027. 6+3335, to 
a model excluding the IFGL source and including two 
sources instead: J2025. 1+3342 and J2028.3+3333. The 
LRT favors the 2-source hypothesis with x^/dof = 282/7, 
corresponding to ~ 16 a. 

All results were cross-checked by an independent anal- 
ysis, and are also found in agreement with a previous 
analy sis of the same fie ld of view including 27 months of 
data ( [Errando fc Ka ra 2010). 

In the following we describe more detailed results for 
these two individual gamma-ray sources, as well as IFGL 
2015.7+3708, all of which are also summarized in Table[2j 
All reported errors are statistical only. The estimated 
systematic uncertainty on t he flux is 5% at 50 MeV, 
increasing to 20% at 10 GeV ( |Abdo et al.]|2010b[ ). 



2.1. IFGL J2015. 7+3708 

In this analysis IFGL J2015.7-h3708 is detected with 
TS — 808, equivalent to approximately 28<j. The lo- 
cation of the gamma-ray source was recalculated us- 
ing 31 months of data to be RA = 20'U5'°33.6% DEC 
= 37° 10' 12" with 95% error circle of 1.'8 that in- 
cludes B2013-f370. The time-averaged spectrum is well 
described with a power-law with photon index F — 
2.57 ± 0.02, and a measured integral flux above 1 GeV of 
(9.26 ± 0.26) X 10~^cm-^s-^. Our result shows a softer 

2.09±0.11, 



spectrum than the on e found in EGRET (F 
Hartman et al.||1999[ ). 



2.2. J2025. 1+3342 

In this analysis of 31 months of Fermi-LAT data we 
resolve a gamma-ray source, J2025. 1+3342, towards the 
direction of IFGL J2027.6+3335. It is located at RA 
= 20'^25'"02.4", DEC = 33°42'00" with 95% error circle 



' All additional sources have been reported in the 2FGL catalog 



http : //fermi .gsfc.nasa. gov/ssc/data/ analysis/user/ 
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10' 



in 10"' 



10' 



J2028.3+3333 



- Exponential cutoff 
Log parabola 

- Power law 




TABLE 3 
Swift-XKT Observation Log 



Energy [GeV] 



Fig. 2. — Gamma-ray spectrum of J2028. 3+3333 . The photon 
distribution The best-fit descriptions of the photon distribution in 
the whole energy range The data was fit to a simple power-law, log- 
parabola and a power-law with exponential cut-off. The spectrum 
is best described by a power-law with exponential cutoff. 

of 3' that includes the blazar B2023+336. J2025.1+3342 
is detected with TS — 289, equivalent to approximately 
17a. The energy spectrum is well described by a power- 
law with a photon index F = 2.94±0.03 and integral flux 
above 1 GeV of (3.25 ± 0.14) x lO-^cm-^s"!. 

2.3. J2028. 3+3333 

In addition to J2025. 1-1-3342 we resolve another, 
brighter gamma-ray emitter, J2028. 3-1-3333, which is 
likely the primary contributor to IFGL J2027. 6-1-3335. 
The second, brighter source, J2028.3-f 3333, is detected 
with TS = 1112, or approxim ately 33cr. As indi- 
cated ^ 



Errando & Kara (20To| the energy spectrum 
is best fit by a power-law with exponential cutoff where 
r = 0.97 ± 0.05, = 1.37 ± 0.03 GeV and integral 
flux above 1 GeV of (10.74 ±0.38) x lO-^cm-^s"!. A 
LRT comparing a power-law, exponential cutoff and log- 
parabola spectral fit favors the exponential cutoff with a 
X^/dof= 159/1 (> 12(t) with respect to a simple power- 
law (Figure [2]) . This source does not exhibit variability 
in gamma-rays. A test on a 7-day binned light curve 
(not shown) gives a probability of the source being vari- 
able of P = 3%. 

2.4. Comparison with Fermi 2FGL catalog 

After our gamma-ray analysis was completed, a first 
version of the Jermi-LAT 2 -year catalog (2FGL) was re- 
leased ( |Abdo et a"L||2011c[ ). The 2FGL catalog, which 
integrates over 24 nionths of LAT data, confirms the re- 
sults obtained for IFGL J2015.7-h3708, J2025. 1+3342, 
and J2028. 3+3333, which are respectively associated 
with 2FGL J2015.6-f3709, 2FGL J2025. 1+3341, and 
2FGL J2028.3+3332. The location and spectral index 
agrees within statistical errors for the three sources. Dif- 
ferences up to 20% in fiux normalization are observed 
in IFGL J2015.7-H3708 and J2025.1+3342, which could 
be enhanced by the different integration time and the 
fact that both sources are variable. Another difference 
is that 2FGL lists a LogParabola function as best de- 
scription for the energy spectrum of the three discussed 
sources, while our analysis preferred power-law spec- 
tra for the two blazars. This difference might be ex- 
plained by the different energy cuts, as it is observed 





Observation ID 


Start time (UT) 


XRT Exposure 






[yyyy-mm-dd hh:mm:ss] 


[sec] 




00035639001 


2006-07-12 05:29:01 


2679.82 




00035639002 


2006-11-12 16:33:01 


4551.96 




00035639003 


2006-11-17 04:09:01 


7449.10 




00035639004 


2006-11-24 04:50:01 


203.65 




00041471001 


2010-08-05 18:18:48 


994.50 




00041471002 


2010-08-06 01:13:38 


7424.83 




00041471003 


2010-08-22 05:35:00 


4360.92 




00041471004 


2010-08-30 01:08:01 


4216.07 



that the biggest deviation from a power-law spectrum 
in the 2FGL sources happens to be in the 100-300 MeV 
bin, which was excluded from our analysis. The other re- 
ported discrepancies could be caused by multiple analysis 
differences: models for the Galactic and isotropic diffuse 
emission, time periods covered (31 months in this study 
vs. 24 in the 2FGL catalog), event calibrations (P6V3 
vs. P7V6), instrument response functions, and energy 
cuts {E > 300 MeV vs. E > 100 MeV), or versions of the 
ScienceTools. Finally, we note that 2FGL J2015.6-H3709 
and 2FGL J2025. 1-1-3341 are respectively associated to 
the radio counterparts of B2013+370 and B2023+336, 
while 2FGL J2028.3+3332 is left unassociated. 

3. X-RAY AND OPTICAL OBSERVATIONS 

The Swift Gamma-Ray Explorer is designed to make 
prompt observations of Gamma-ray Bursts and after- 
glows, but also performs pointed observations. For this 
analysis, we use data from the Swift X-ray T elescope 
(XRT ) , which is sensitive from 0.2 to 10 keV (Gehrels 
e raLl| 2004). 

During August 2010, four Swift pointed observations 
of B2013-f370 (IFGL J2015.7+3708) were triggered by 
a Fermi detected flare on 2010 August 4\ Archival data 



from observations taken in 2006 were also analyzed. We 
performed two separate analysis on the 2010 and 2006 
observations. Table [3] summarizes the Swift-XRT obser- 
vations of B2013+370. No Swift data on B2023+336 are 
available. 

The XRT data were processed with standard proce- 
dures in xrtpipeline, adopting the standard filtering 
and screening criteria, using FTOOLS in the Heasoft 
package (v6.10Q For both the 2006 and 2010 sets, we 
summed the data from the four individual observations. 
Due to the low count rate of this weak source (6.0 x 10^^ 
counts s^i in 2010, and 2.3 x 10"^ counts s^^ in 2006), 
we consider only photon co unting data, and furt her se- 
lect XRT event grades 0-12 (|Burrows et al.||2005|). Pile- 
up correction was not rcquirecTTTrom the 2010 and 2006 
counts maps, we extracted the source events from a cir- 
cular region with radius 20 pixels (1 pixel = 2.36"), and 
also extracted background events within a 55-pixel radius 
circle that was offset from the source. 

Ancillary response files were generated with xrtmkarf , 
accounting for the different extraction regions, vignetting 
and PSF corrections. We used current response ma- 
trix file vOll available from HEASARC Calibration 

'http : / /f ermisky .blogspot . com/ 
http : //heasarc . gsf c.nasa. gov/^ 
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Databasep^ The spectra were rebinned in order to have 
at least 20 counts per energy bin, and the spectral fitting 
were completed using XSPEC vl2.6.0q. 

The X-ray spectra were fit by a power-law F{E) = 
KE^^ , absorbed by the Galactic column density, A^^^' ~ 



x10° 



1.7 X 10^^ cm- 



which was determined from the 



data. The spectrum from the August 2010 obser- 
vations resulted in a photon index T = 1.77^q'20' 
and normalization at 1 keV oi K ^ l-46lQ'3g ^ 
10^'^ keV^^ cm^^ s~^. The spectrum from the 2006 ob- 
servations showed a similar photon index and higher flux 
(r = 1.23i°i^, K = 3.26i;jJ X 10-3 keV'^ cm'^ s~^). 

Archival optical observations, imaged on 2004 July 
10, were available for the position of B2023-f-336 from 
the 2.4m Hiltner Telescope of the MDM Observatory 
with a thinned, backside-illuminated SITc CCD having 
a spatial scale of 0'.'275 per 24/x pixel. Three 5 minute 
exposures in each of the B, R, and / filters were ob- 
tained in seeing of O'.'S, and were combined. An as- 
trometri c solution for the ima ges was derived from the 
UCAC2 ( Zacharias et al.|2004 ). The optical counterpart 
of B2023-I-336 was clearly detected in each band, at po- 
sition R.A.=20'^25'"10."85, decl.=-h33°43'00"3 (J2000.0), 
which is consistent with the radio p osition. Magnitudes 
were calibrated with Landolt ( 1992 ) standard stars. In- 
cluding both statistical and systematic errors, we find 
B = 24.12±0.14,i? = 21.07±0.04, and / = 19.43±0.04. 

Galactic extinction is a significant factor in the direc- 
tion of B2023+336, = (73°. 13, -2°. 37 ). The Hy- 



drogen c olumn density (iVn = 6.27 x 10 Kalberla et 



al.|2005 1 was converted to optical extinction,^ 



using the linear relation described in Giiver & Oumlzel 
( 2009 1, and scaled to the selected bandpasses using ratios 
given in |Schlegel et al., ( ,1998^ . The extinction-corrected 
optical fluxes tor the B, K and / band can be found in 
Table H 

4. RADIO OBSERVATIONS 

Contemporaneous 15 GHz radio light curves were ob- 
tained at the radio positions of the blazars B2013+370 
and B2023-I-336 with the 40-m telescope at the Owens 
Valley Radio Observatory (OVRO). Since late 2007, 
about a year before the launch of Fermi, a large-scale, 
fast-cadence 15 GHz radio monito ring program has been 
ongoing using the 40-m telescope ( Richards et al.pOll ). 
The OVRO 40-m program began with the 1158 nortfi- 
ern {5 > —20°) sources from the Ca ndidate Gamma-ray 

and now includes 



2008), 



Blazar Survey (Healey et al. 
over 1500 sources in cluding ail the nort hern Fermi ILAC 
blazar associations (Abdo et al.|2010c|. The two sources 



in this analysis were added to the OVRO 40-m program 
in August, 2008. Program sources are each observed 
twice per week in total intensity with about 4 mJy (min- 
imum) and 3% (typical) un certainty. The abs olute flux 
density scale is based on the Baars et al. ( 19771 value for 
3C 286, 3.44 Jy at 15 GHz, with about K^Tscale uncer- 
tainty which is not reflected in the error bars presented 
here. 

Radio light curves of B2013+370 and B2023-I-336 are 
respectively shown in Figure [3] and [4j and the time- 
averaged radio fluxes are reported in Table [5] and |4] 



0.3 



0.2 



StL 0.1 



1FGL J201 5.7+3708 {E>300 MeV) 




B201 3+370 (15 GHz) 



54600 



54800 



55000 55200 
time [MJD] 



55400 



55600 



Fig. 3. — Top: Gamma-ray light curve at i? > 300 MeV in 7-day 
bins for IFGL J2015. 7+3708 (B2013-I-370). Open circles indicate 
fluxes for time periods where the source was detected with TS < 
4, and grey triangles represent 90% c.l. upper limits (TS < 1). 
Bottom: OVRO radio light curve at 15 GHz for same 31 month 
period. 

5. LOCATION OF THE GAMMA-RAY SOURCE 

A detailed analysis of the Fermi-LAT data in the di- 
rection of B2013+370 and B2023-I-336 was presented in 
Section [2j The resolved gamma-ray sources and their lo- 
cations with respect to the radio positions of the blazars 
are shown in Figure [T] 

The possible association of the EGRET unidentified 
source 3EG J2016-h3657 with the blazar B20 1 3-K37 was 
extensively discussed 



(2000) and 



Mukherjee et al. 
Halpern et al.] |2001|). The authors examined X-ray im- 
ages of the EGRET error box and conducted optical 
spectroscopy of all candidate counterparts, finding the 
blazar B2013-I-370 the only probable counterpart. With 
an 11-month exposure, the Fermi-hAT collaborati on re- 
ported a gam ma-ray excess (IFGL J2015. 7+3708, lAbdo 



et al. 



associated with 3EG J2016+3657, 
-370 lies just outside the IFGL 95% error 



2010^ 

ElaiaJ'B^mS 

ellipse. A slight displacement of the Fermi reconstructed 
position could be caused by the high level Galactic dif- 
fuse emission in the region and the presence of a bright 
gamma-ray pulsar o nly 1.2° away (PSR J2021-)-3651, 
Abdo et al.| |2010d|). In the analysis of Fermi-LAT 



http : //heasELTC ■ nasa . gov/ docs/heasarc/caldb| 



data presented here, we recalculated the location of the 
gamma-ray source with 3 times better photon statistics 
(see Tablej2]), and found the radio position of B2013-I-370 
compatible with the reconstructed gamma-ray location 
(Figure [T]). The newly available Fermi 2-year catalog 
confirms this association. Recently, the VERITAS col- 
laboration reported the detection of a VHE gamma-ray 
source compatible with th e nearby p ulsar wind nebula 
CTB 87 at > 650 GeV ( |Alm|[20Tl| ). The emission is 
steady and significantly displaced from B20 13-1-370, sug- 
gesting a different origin. 
The blazar B2023-I-336 was proposed as the 
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TABLE 4 
B2023+336 Spectral Points 



Tn G'f' Tvi d 








[Hzl 

L J 


[10~^2 erg s"'^ cm~2] 


OVRO 


1.5x10^" 


0.428 ±0.013 


MDM 


3.75x10^'' 


0.671 




4.61x10" 


0.420 




6.82x10" 


0.229 


_Ferm j-LAT 


1.02x10^^ 


21.2 ± 2.5 




3.25x10^3 


6.19 ± 1.13 




1.04x1024 


3.31±0.79 




3.32x1024 


< 1.01 




1.06x10^5 


< 1.29 



xlO" 



E 

^ 0.10 



J202S.U3342 (E>300 MeV) 




t / 



54900 54950 



62023+336(15 GHz) 



A. ' • 



54600 54800 55000 55200 55400 55600 

time [MJD] 

Fig. 4. — Top: Gamma-ray light curve at i? > 300 MeV in 14-day 
bins for J2025.H-3342 (B2023-I-336). Open circles indicate fluxes 
for time periods where the source was detected with TS < 4, and 
grey triangles represent 90% c.l. upper limits {TS < 1). Bottom: 
OVRO radio light curve at 15 GHz for same 31 month period. The 
inset shows a zoom-in of the radio light curve around MJD 54930, 
when a flux modulation tentatively associated with an extreme 
scattering event is observed. 

counterpart for the EGRET unidentified source 



3EG J2027 -H3429 ( [Sowards-Emmerd et al.||2003[ [Sguera 
et al. 2004 1 . The Fermi-LAL' ll-month catalog reports 
an unidentified source (IFGL J2027. 6-1-3335) outside the 
EGRET error eUipse. The position of the IFGL source 
is 32' away from the radio position of B2023-f336. In our 
analysis presented in Section [2] we resolve a component 
of IFGL J2027.6-H3335: J2025.1-I-3342, spatially associ- 
ated with the radio location of B2023-f336 (Figure [T]). 
The recent release of the Fermi 2-year catalog confirms 
this association. 

6. FLUX VARIABILITY AND CORRELATIONS 

Variability in gamma-ray sources allows an identifica- 
tion of the correct counterpart when correlated multi- 
wavelength variability is found. In addition, stand-alone 
gamma-ray variability can also support the affiliation of 
an unidentified source to an intrinsically variable source 



class (see, e.g. Nolan et al. 2003). Flux variability is 
one of the main observati onal characteristics of blazars 
in the gamma-ray band (Abdo et al. 2Q10f), whereas 
only a handful of non-AGW sources have shown s ignifi- 
cant g amma-r ay variability: 5 gam r na-ray binar ies (Abdo] 



et al. 



2010b[ [Corbet et al.ll20TT) [Abdo et al ,2011bj f 
Jygni (|S)do ct al. 2010c), and the Grab'iNebulf 



( [Tavani et alpOll, ,Abdo et al.||2011a ). 

I'he light curves of the gamrna-ray sources associated 
with B2013-h370 and B2023-h336 exhibit some degree of 
variability. We estimated the gamma-ray variability with 
a test on th e light curve, i ncludin g all flux points with 
TS* > 1 as in ^ 
light curve of 



Abdo et al. (2009b). The 7-day binned 

muimFT+ms (B2013+370) shown 

in Figure [3] exhibits clear flaring episodes, and signifi- 
cant evidence of the flux being variable (P > 99.9%). 
J2025. 1-^3342 (B2023+336) is a weaker gamma-ray emit- 
ter. We calculate a probability of its flux being variable of 
P = 99% based on the 14-day binned light curve shown 
in Figure 4] This cannot be considered conclusive evi- 
dence of the source being variable. However, as noted in 
Abdo et al. (2010c), given the same fractional flux vari- 



ation it is harder to find significant variability in faint 
sources like J2025. 1-1-3342 than it is for bright sources 
because of the strong dependence of the variability test 
on the statistical flux uncertainties. 

A possible correlation between the gamma-ray and ra- 
dio light curves is studie d by constructing a discr ete cor- 
relation function (DCF, Edelson & Krolik 1988) shown 
in Figure [Sj The statistical significance of the cross- 
correlation is investigated with Monte Carlo simulations. 
It is assumed that the light curves can be described as 
noise processes with a power-law power spectral density 
(oc /^). The adopted power-la w exponents are - 1.5 for 
gamma -rays jAbd o et al.|2010f| and -2.0 for radio ( IChaF 
terjee et al. 2008). Using the recipe prescribed by Tim 



Koenig ^ ( .19951 , 
th a dai 



mer 

is generated witTi 



a pair of independent light curves 
ily sample rate. These light curves 



are sampled at the same times as the actual radio and 
gamma-ray light curves and the values are perturbed by 
adding Gaussian white noise with zero mean and vari- 
ance that matches the observational errors. The result- 
ing mock light curve s are then cross-correlate d using the 
method described in Edelson & Krolik (1988). The pro- 
cess is repeated 10000 times, and the results provide an 
estimate of the distribution of random cross-correlations 
for each time lag. 

The results of the cross-correlation tests are summa- 
rized in Figure [5] which shows the DCF for the data 
as black points and the 1-cr (red), 2-a (yellow) and 3-cr 
(green) contours for the distribution of random cross- 
correlations. Using these values we find no highly signif- 
icant cross-correlation. The most prominent peaks are 
located at -10 days (radio lagging) with an 88.0% sig- 
nificance for B2013-I-370 and a peak at -240 days with 
a 96.7% significance for B2023-I-336. Additionally, the 
auto-correlation functions (not shown) were calculated 
separately for the gamma-ray and radio light curves, 
without finding any significant peak that could suggest 
periodic behavior. 

Although a statistically significant correlation between 
the gamma-ray and radio light curves would provide a 
definite identification, such identiflcations based on cor- 
related variability have only been established for very 
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B2013 + 370 B2023 + 336 




T[days] T[days] 



Fig. 5. — Discrete correlation function (DCF) between the gamma-ray and radio light curves for B2013-I-370 (left) and B2023+336 (right). 
Positive time lags indicate that radio activity precedes activity in gamma-rays. The significance of the correlations is illustrated by l-cr 
red, 2-(T yellow and 3-(t green dotted lines calculated from simulated uncorrelated data assuming a simple power-law power spectral density 
with exponent -2.0 for radio and -1.5 for gamma-ray light curves. The inset DCF for B2013+370 during the gamma-ray flare in MJD 
54750-55000. 



TABLE 5 
B2013+370 Spectral Points 







Average 




2010 August 


Instrument 


Frequency 




Frequency 






[Hz] 


[10^1^ erg s^i cm^^] 


[Hz] 


[10^12 erg s~i cm"2] 


OVRO 


1.5x10^0 


0.485 ±0.01 


1.5x10^" 


0.453 ±0.009 


Swift-XRT 


5.11x1017 


1.43 ±0.10 


5.84x101'^ 


0.74 ± 0.09 




8.42x10^7 


2.16 ±0.15 


7.10x101'^ 


2.15 ±0.14 




1.19x10^** 


3.05 ±0.22 


9.66x101'^ 


3.05 ±0.24 




1.51x101** 


3.72 ±0.27 






Fermi-LAT 


9.81x10^^ 


36.1 ± 10.3 


1.05x10^3 


68.8 ± 13.3 




2.25x10^3 


29. 1± 3.6 


3.35x10^3 


47.7 ± 10.4 




5.16x10^3 


15.8 ±2.1 


1.07x10^4 


24.0 ±9.7 




1.18x10^3 


7.87 ± 1.54 


3.42x10^4 


< 10.9 




2.71x10^* 


3.64 ± 1.29 


1.09xl02"' 


< 33.2 




6.22xl02'i 


< 4.53 








1.43x10^5 


< 8.35 







bright gamma-ray blazars during short-lived flares (see, 



e.g 


. IWehrle 


et al. 


111998 


al. 


2010a). 


If we 


only i 



correlation peak at -10 days lag clearly increases, as 
shown in Figure [Sj Assessing the statistical significance 
for very short light curves is difficult because the noise 
properties used for the statistical test are derived from 
longer time series and might not be appropriate for short 
periods of time, especially when these have been selected 
because of unusual source activity, a prominent gamma- 
ray flare in this case. The fact that the significance de- 
creases as more data are included indicates that no ro- 
bust claims about the physical significance of the appar- 
ent correlations can be made using time series dominated 
by a single event. 

Intrinsic long-term flux correlations might be washed 
out by the presence of extrin sic variability affecting only 



one of the frequency bands. Spangler & Cordes ( 1998 1 
showed that observations of extragalactic sources m the 
line of sight to B2013+370 and near B20234-336 show 
heavy scattering due to interstellar plasma turbulence 



in the direction of the Cygnus OBI association. This 
mechanism would produce extrinsic variability in the ra- 
dio band but not in gamma-rays. Disentangling intrinsic 
and extrinsic variability in the radio light curves is a dif- 
ficult problem, but a particular feature on the light curve 
of B2023-H336 around MJD 54930 strongl y suggests the 
possi bility of an extreme scattering event (Fiedler et al. 
19871. Although extreme scattering everits have only 



been detected at lower radio frequencies, the location 
of the source at low Galactic latitude and the charac- 
teristic symmetric flux decrease expected when an inter- 
stellar plasma lens moves through the line of sight, with 
the small flux increases at both ends of the event being 
caused by focusing of radio waves at the edg es of the 
plasma cloud, strongly support this hypothesis (Clegg et 
,al.^^l998) . To explore the effect of excluding the candi- 
date extreme scattering event from the radio light curve 
of B2023+336, we use the same parameters as for the 
original analysis but exclude the radio data from MJD 
54900 through MJD 54975. We find a slightly higher sig- 
nificance for the cross-correlation of 98.8% at the same 
time lag. This supports the suggestion that radio-specific 



8 



log(E/eV) log(E/eV) 

-6 -4 -2 2 4 6 8 10 12 -6 -4 -2 2 4 6 8 10 12 




8 10 12 14 16 18 20 22 24 26 8 10 12 14 16 18 20 22 24 26 

log(v/Hz) log(v/Hz) 



Fig. 6.— Spectral energy distribution for B2013+370. Left: Filled points (this work) correspond to OVRO (15 GHz), Swift X RT (2-10 
keV) and F ermi-LAT data {E > 300 MeV). Open circles show archival radio data obtained from NED, R-band magnitudes from 'Halpernj 
|et al.| ([2001), and an _E > 350 GeV upper limit from Fcgan et al. (2005). Right: The blue (average) points consist of XRT data from 2U0S7 
and an average over the 31 months of Fermi and OVKU data. Ihe red (August 2010) points consist of the Fermi and OVRO data from the 
entire month of August 2010, when there was a flare in gamma-rays, which triggered four XRT observations. Solid colored lines indicate 
the total predicted SED by the different tested models assuming a redshift of z = 0.1 (left) and z = 0.2 (right). Dashed lines show the 
Synchrotron emission, dotted lines the SSC component, and dot-dashed lines the EC. 



TABLE 6 

SED MODELING PARAMETERS 



Parameter 


Symbol 


SSC 


SSC+ECi 


SSC-I-EC2 


Electron power 


Le [erg s~i] 


4.6 X 10'*'' 


4.7 X 10'*"' 


1.0 X lO"*" 


Electron low-energy cutoff 


'ymin 


4.0 X 10^ 


1.5 X 10^ 


1.4 X 10^ 


Electron high-energy cutoff 


^max 


3.0 X 10^ 


3.0 X 10^ 


2.5 X 10"* 


Electron injection index 


Qe 


3.2 


3.1 


2.8 


Equipartition fraction 


Lb/Lb 


1.58 X 10-5 


6.49 X 10-2 


1.00 


Blob radius 


Rb [cm] 


7.0 X W^"^ 


6.0 X 10^'' 


1.1 X 10^® 


Magnetic field 


B [G] 


0.001 


0.1 


1.0 


Bulk Lorentz factor 


F 


20 


15 


15 


Observing jet angle 


dobs [deg] 


2.87 


3.82 


3.82 


Redshift (assumed) 


z 


0.1 


0.1 


0.1 


External radiation energy density 


Uext [erg cm"^] 




2.0 X lO-'' 


2.5 X 10-^ 


External radiation temperature 


Text [K] 




10^ 


10^ 



extrinsic effects reduce the observable cross-correlation. 

7. SPECTRAL ENERGY DISTRIBUTION 

The SED for B2013-F370 (Figure [6) data points in 
Table l5| shows a two component structure characteris- 
tic of blazars. The relatively hard X-ray spectral index 
{ax = 1-2 — 1.8) and the dominance of the gamma-ray 
power output are typically observed in low-frequency- 
peaked B L Lacs (LBLs ) or in flat spectrum radio quasars 
(FSRQs, [Fossati et alT]|l9 98). 

We attempted to describe the multi-wavelength SED 
of B2013-I-3 70 with the one-zone lepton ic e mission model 
described in[B6ttchcr & Chiang ("2002) and Acciari et al. 
( |20()9| ). The model-predicted SEDs arc shown in i'lg- 
ure |m and the assumed model parameters are summa- 
rizea in Table |6] The redshift of the source was assumed 
to be 0.1 throughout the modeling proce ss. A pure syn 



of the SED and requires a very low magnetic field, result- 
ing in a strongly particle dominated jet [Lb/ L^. ^ 10^^, 
with Lb and L^ being the magnetic and particle power). 
This is not surprising: SSC models typically fail to repro- 
duce the high-energy component of Compton dominated 
blazars (LBLs, and specially FSRQs) unless the mag- 
netic field is reduced to unphysically low values. This 
is addressed in external C ompton mo dels (EC, e.g. |Der- 



chrotron s elf-Comp ton model (SSC, e.g. IVIaraschi et al. 
[1992^ Bloom fc ]VIarschcr„1996| gives a poor description 



mer k Schhckeiser,,1993| [Sikora fc Madcjski|[2000[ r^ 
adding a population ot low energy photons external to 
the jet. FSRQs are often described with leptonic models 
where SSC is dominating in the X -ray band and EC ex - 
plains the gamma-ray output (e.g. Hartman et al.|2001 ). 
We explore this scenario by adding an external radiation 
field with temperature T^xt — 1000 K. The resulting set 
of parameters (SSC-l-ECi) is closer to the typical val- 
ues encountered in blazar SED modeling, but still needs 
low sub-equipartition magnetic fields and fails to repro- 
duce the observed X-ray spectral slope. The best fit 
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Fig. 7. — Spectral energy distribution for B2023+336 using 
OVRO (15 GHz), MDM (B, R, and I bands), and Fermi data 
{E > 300 MeV). Open points indicate archival data from NED, 
and from two observ ations taken with BeppoSAX (1-100 kcV) from 
[Sguera et ar| | |2004[ l. 

(SSC+EC2) is achieved with both X-rays and gamma- 
rays being dominated by EC emission. Parameters in 
exact equipartition {Lg/Lf, = 1) can be achieved in this 
case, requiring a rather low temperature of the external 
radiation field {Text = fOOK) that could originate from 
cold dust. We note that the relatively low frequencies of 
the photons produced in the first-order SSC (~ 10^^ Hz) 
make the second-order SSC (peaking at ^ 10^^ Hz) be- 
come efhcient. The lack of an £' > 50GeV detection, 
where the gamma-ray absorption by extragalactic back- 
ground light is strongly redshift dependent, makes the 
assumed z = 0.1 a non-critical model parameter. Sim- 
ilar best-fit solutions close to equipartition were found 
assuming 2; up to 1.5. A similar set of model parame- 
ters (assuming z — 0.2 this time) can also describe the 
average and 2010 August states shown in Figure |6] 

The SED for B2023-H336, shown in Figure[7l shows the 
same blazar-like features discussed for B2013-I-370: two- 
component structure, hard X-ray spectrum, and dom- 
inant gamma-ray power. However, only archival non- 
contemporaneous optical and X-ray data are available. 
Given that B2023-f336 is known to be variable at all 
frequencies, we don not attempt to model this source 
because of the lack of a simultaneous SED. 

8. DISCUSSION AND CONCLUSIONS 

Spatial association and the observed variability in 
the gamma-ray and radio bands allow us to establish 
a firm association between B2013-I-370 and the previ- 
ously unidentified gamma-ray source 3EG J2016+3657 
(IFGL J2 15.7-f- 3708), con firming the t hesis of iMukher 



jee et al. (2000) and Halpern et al. (2001). A com 



piled SED, aHffing newly analyzed OVEnT Swift-XRT 
and Fermi-LAT data shows a two component structure 
that further supports the blazar association. The SED 
was successfully modeled with a leptonic one-zone model 
with both X-ray and gamma-ray power dominated by 



external Comptonization of a low-temperature external 
radiation field. The gamma-ray dominated SED, hard 
X-ray spectrum, and preference for EC models point to- 
wards B2013-h370 being an LBL or an FSRQ. 

We also find strong evidence of gamma-ray emission 
from the blazar B2023+336, that was likely contributing 
to the unidentified gamma-ray source 3EG J2027+3429 
IFGL J2027.6+3335 ), as hypothesized by [Sowards" 



Emmerd et"ar| ( |2003| and [Sguera et aT] ( |20(J4p . Uur 

analysis of the available Fermi-LAi' data resolves two 



analysis of the available Fermi-LAL data resolves two 
independent gamma-ray sources of very different charac- 
ter. J2025. 1+3342 is spatially associated with the blazar 
B2023-I-336 and shows a hint of variability. The tentative 
observation of an extreme scattering event in the radio 
light curve of B2023-f 336 further supports the hypothesis 
of it being extragalactic in origin, the flux decrease being 
caused by a cloud of plasma in the galactic interstellar 
medium crossing the line of sight. 

The other resolved component of 3EG J2027+3429 
(J2028. 3-1-3333) is characterized as a steady gamma-ray 
emitter with an exponentially cutoff energy spectrum, 
qualities which are typical of LAT detecte d pulsars. This 
result has recently been confirmed by [Pletsch et al.| 
(20111, who report the discovery of nine previously un- 



known gamma-ray pulsars using a blind s earch method. 
This d iscov ery supports the hy pothesis of Romero et al. 
( |1999| and [Zhang et al] ( [200b( ) that a pulsar could be 
responsible tor the gamnia-ray excess 3EG J2027+3429, 
or at least part of it. 

Optical spectroscopy of the studied objects is compli- 
cated due to their low optical luminosity and the ex- 
treme dust absorption in the optical band characteristic 
of low Galactic latitudes. Future observations in large- 
aperture optical telescopes could establish a definitive 
spectroscopical identification of both blazars. 

This work has been made possible for the first time 
because of the continuous, well sampled light curves in 
the gamma-ray and radio bands obtained over 2.5 years 
by Fermi and the fast-cadence OVRO blazar monitoring 
program. In completing this analysis, we establish firm 
associations for the blazars B2013+370 and B2023+336 
with previously unidentified EGRET and Fermi-hAT 
sources. This confirms the gamma-ray band as likely the 
best tool to identify blazars located behind the Galactic 
plane, which are heavily absorbed at other wavelengths, 
and supports the hypothesis that a number unidentified 
gamma-ray sources at low Galactic latitudes are indeed 
of extragalactic origin. 
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